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MHD  SIMULATION  OF  A  COMET  MAGNETOSPHERE 

Introduction 

The  approaching  1986  apparition  of  Comet  P/Ealley  has  generated  a 
great  deal  of  scientific  interest  and  inquiry  concerning  the  interaction  of 
the  solar  wind  with  cometary  plasmas.  tfith  the  planned  space  probe 
missions  to  make  direct  measurements  of  Comet  P/Halley  plasma  properties, 
there  is  considerable  Interest  to  make  predictions  concerning  the  solar 
iflnd-comet  interaction.  This  paper  presents  the  results  of  a  numerical 
simulation  of  this  interaction  using  magnetohydrodynamlc  (MED)  theory.  Our 
results  can  be  described  in  terms  of  a  cometary  magnetosphere  owing  to 
similarities  between  the  comet  and  the  earth's  magnetosphere;  particularly 
a  strong  magnetic  field  separating  the  cometary  body  and  the  solar  wind  on 
the  bow,  and  a  very  long  tenuous  magnetic  tail  with  oppositely  directed 
field  in  two  lobes  separated  by  a  plasma  sheet  and  neutral  field  region. 
However,  the  analogy  between  the  earth  and  a  comet  is  limited  as  the  plasma 
dynamical  properties  are  considerably  different. 

The  original  concept  of  a  cometary  magnetosphere  was  developed  by 
Blermann  (1951),  who  showed  that  the  main  features  of  comet  plasma  tails 
could  be  determined  by  Interaction  with  an  as  yet  poorly  understood  solar 
wind.  He  also  demonstrated  that  comets  could  serve  a  useful  purpose  as 
probes  of  the  solar  wind.  The  next  major  step  in  understanding  the  solar 
wind  cometary  interaction  came  from  Alfven  (1957)  who  discussed  the  effect 
of  the  Interplanetary  magnetic  field  (IMF)  on  the  interaction.  Alfven 
pointed  out  that  the  IMF  would  be  "hung  up”  in  the  head  of  the  comet  and 
the  solar  wind  would  drag  the  field  into  a  long  folding  tall,  quite  like 
the  folding  of  an  umbrella.  Subsequent  research  by  Marochnik  ( 1963a, b)  and 
Axford  (1964)  showed  that  the  solar  wind  interaction  with  comets  could  be 
described  by  MHD  theory. 

Many  continuing  efforts  considered  details  of  the  solar  wind  comet 
Interation  in  the  bow  region,  particularly  Blermann  et  al.  (1967),  Ioffe 
(1966a, b;  1968a, b),  Wallis  (1973),  Brosowskl  and  Wegmann  (1973),  and  Wallis 
and  Dryer  (1976).  These  works  treated  the  existence  of  a  shock  in  the 
solar  wind  preceedlng  the  comet  and  a  contact  discontinuity  enclosing 
plasma  of  purely  cometary  origin.  Much  of  this  work  and  a  summary  of 
results  is  discussed  in  excellent  reviews  by  Mendls  and  Ip  (1977),  Brandt 
and  Mendls  (1979),  Breus  (1982),  Mendls  and  Houpls  (1982)  and  Ip  and  Axford 
MasMCftpt  approvsd  January  31, 19M. 


(1982).  An  excellent  review  of  the  theory  of  Instabilities  of  cometary  in 
tails  has  been  presented  in  Ershkovlch  (1980).  In  a  series  of  recent 
papers  by  Niedner  and  Brandt  (1978,  1979),  Neidner  et.  al.  (1981),  and 
Brandt  (1982),  the  interaction  of  a  dynamically  changing  solar  wind  with 
comet  magnetospheres  has  been  studied,  particularly  concerning  comet  ion 
tall  dynamics  and  tall  disconnection  events. 

Recently,  numerical  simulations  of  the  solar  wind  comet  interaction 
have  been  presented  by  Schmidt  and  Wegmann  (1980,  1981).  In  these  works 
they  studied  plasma  flow  between  a  hydrodynamic  shock  at  some  distance  from 
the  comet  and  an  arbitrary  fixed  contact  surface  near  the  comet  nucleus. 
Their  field,  flow  and  plasma  density  structure  are  appropriate  for  the 
inner  regions  of  a  comet  magnetosphere.  Unfortunately  their  numerical 
algorithm  was  unable  to  maintain  divergence  B^  identically  zero  over 
substantial  regions  of  the  simulation.  Brackbill  and  Barnes  (1980) 

demonstrated  the  serious  difficulties  in  numerical  simulations,  not  only  in 
magnetic  topology  but  also  in  plasma  diagnostics,  which  can  arise  caused  by 
a  divergent  magnetic  field. 

In  this  paper  we  study  the  solar  wind  interaction  in  the  outer  regions 
of  a  cometary  magnetosphere.  The  comet  acts  as  a  source  for  production  of 
plasma  in  the  solar  wind.  It  is  seen  that,  as  Alfven  proposed,  the  IMF 
does  "hang  up"  in  the  region  of  the  comet.  The  magnetic  field  is  pulled 
into  a  long  tall  confining  the  cometary  plasma.  We  demonstrate  the 
configuration  of  field  and  plasma  for  a  steady  solar  wind  and  IMF.  We  also 
'  explore  the  effect  of  rotational  discontinuities  in  the  IMF  on  the  comet 

magnetosphere. 

Mathematical  and  Numerical  Model 

We  consider  a  solution  of  the  ideal  MHD  equations  on  a  3  dimensional 
cartesian  mesh  (x,y,z).  The  equations  solved  numerically  have  the  form  (in 
gausslan  units): 
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In  these  equations,  p  Is  the  mass  density,  ^  is  the  plasma  flow  velocity,  P 
Is  the  plasma  pressure,  ^  Is  the  electric  Intensity,  JB  is  the  magnetic 
induction,  is  the  electric  current  density,  c  is  the  speed  of  light,  y  is 
the  ratio  of  specific  heats,  A  is  the  comet  plasma  production  rate,  and  Q 
is  an  artificial  viscosity  which  is  used  to  numerically  conserve  energy  in 
the  plasma  (for  example  Q  is  formulated  to  meet  jump  conditions  at  shocks). 

The  method  of  solution  is  as  follows:  .  Equations  (1),  (2),  (3),  and 
(5)  are  explicitly  integrated  in  time.  Equation  (4)  and  (6)  are  used  to 
substitute  for  dependent  variables.  The  numerical  algorithm  used  is  the 
partial  donor  cell  method,  PDM,  (Hain,1978).  It  is  a  nonlinear  hybrid 
algorithm  of  the  shock  capturing  t3rpe.  The  PDM  algorithm  has  been  run  on 
numerous  test  problems  with  excellent  results.  By  careful  treatment  of 
equation  (5),  divergence  B_  is  identically  zero  throughout  the  simulation, 
thereby  eliminating  possible  serious  numerical  error  (Brackblll  and  Barnes, 
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1980).  The  numerical  method  has  been  used  previously  to  simulate  the 
earth's  magnetosphere  Brecht  et  al.  (1981,  1982)  and  has  shown  results 

which  agree  well  with  theory  and  data. 


For  the  comet  simulation  we  consider  an  Initial  value  problem  and 
solve  the  equations  In  time  to  a  quasi  steady  condition.  The  simulation  Is 
Initialized  on  a  cartesian  mesh  which  extends  from  1.4  x  10^  km  In  front  of 
the  comet  to  5  x  10^  km  In  the  antlsolar  direction,  and  to  1.25  x  10^  km  on 
each  side.  A  steady  solar  wind  flows  In  the  x-dlrectlon  having  a  velocity 
of  400  km/ sec,  a  density  of  5  cm  ,  a  temperature  of  5  ev  and  an 
Interplanetary  magnetic  field  of  5y  In  the  z  direction.  These  values  are 
appropriate  to  one  astronomical  unit  distance  from  the  sun.  We  take  the 
cometary  plasma  production  rate  of  Schmidt  and  Wegmann  (1980)  which  Is 
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to  gradually  build  up  the  comet.  The  cometary  plasma  production  parameters 
In  (7)  are:  ■  30  proton  masses,  G  ■  10^®  cm""^,  w  ■  I  km/sec,  <j  ■  10”® 
sec”^,  and  the  radius  from  the  comet,  R  ■  (x^  +  y2  +  In  order  to 
study  the  result  of  asymmetric  plasma  production  due  to  extinction  of  solar 
Ionizing  radiation  we  have  also  multiplied  A  by  the  quantity  (1  -  x/R). 
This  particular  function  was  not  chosen  to  model  a  specific  extinction 
process  but  was  used  to  test  If  extinction  was  Important.  In  either  case, 
the  simulation  settles  down  to  a  steady  state  after  about  12  hours  of  real 
time. 


Results 

Results  have  been  obtained  for  the  structure  of  a  cometary 
magnetosphere  with  both  spherically  symmetric  and  as3nnmetrlc  plasma 
production.  All  results  shown  here  are  for  the  asymetrlc  case  since  the 
two  cases  where  virtually  Indistinguishable,  both  qualitatively  and 
quantitatively.  For  the  simulation  the  solar  wind  flows  In  the  positive  x 
direction.  The  results  show  a  very  distinct  plasma  tall  structure  for  a 
steady  solar  wind  and  IMF.  Results  for  a  rotational  discontinuity  In  the 
IMF  were  also  obtained.  These  results  show  an  Indication  of  a  tall 
disconnection  event  and  also  ray-like  structure  formation. 


We  begin  Che  discussion  of  our  results  by  considering  the  steady  solar 
wind  and  IKF  case.  Figure  1  shows  the  magnetic  field  structure  in  the 

comet  for  the  plane  defined  by  Che  sun-comet  line  and  the  direction  of  Che 
IMF.  Figure  la  shows  contours  of  magnetic  field  strength  and  lb  shows  unit 
vectors  of  the  field  direction.  The  maximum  magnetic  field  strength 
of  48y  occurs  at  a  distance  1  x  10^  km  in  front  of  nucleus.  The  field 
rises  quickly  to  its  maximum  value  in  the  head  and  decays  very  gradually 
into  the  tall  where  at  4  x  10^  km  behind  the  nucleus  the  strength  is  lly. 
The  cometary  plasma  sheet  in  the  tail  has  a  magnetic  field  which  is 
considerably  smaller  than  Che  IMF  value  and  is  very  narrow  (>'  10^  km  in  the 
direction  of  Che  IMF).  The  field  vectors  demonstrate  Che  draping  of  Che 
IMF  over  the  dense  plasma  in  the  cometary  head. 

Figures  2a  and  2b  show  the  flow  velocity  and  the  plasma  pressure, 
respectively,  in  the  plane  of  the  IMF.  The  length  of  the  velocity  vectors 
indicates  flow  speed.  The  most  striking  feature  of  the  flow  pattern  is 
Chat  Che  direction  of  the  flow  is  virtually  unchanged  from  the  solar  wind 
direction  throughout  the  cometary  plasma.  There  is  only  a  very  slight 
deviation  of  flow  which  occurs  at  a  weak  shock  located  2  x  10^  km  in  front 
of  the  nucleus.  The  solar  wind  flow  decelerates  rapidly  in  front  of  and  at 
Che  shock  front  and  then  reaccelerates  very  gradually  behind  the  comet.  In 
the  near  comet  central  tail  the  flow  velocity  is  very  low,  only  about  10 
km/sec.  In  Figure  2b  the  perspective  plot  of  plasma  pressure  shows  a 

number  of  interesting  features.  We  see  the  compresslonal  heating  of  the 
solar  wind  flow  to  a  few  hundred  eV  cm~^  slightly  ahead  of  and  at  the  shock 
front  and  a  cooling  of  the  plasma  behind  the  shock  as  the  flow 
reaccelerates.  There  is  also  a  very  thin  plasma  sheet  which  sits  in  Che 
magnetic  field  neutral  sheet.  In  the  plasma  sheet  the  Increased  pressure 
is  due  primarily  to  large  plasma  densities  and  not  to  high  plasma 

temperatures.  The  one  point  maximum  in  the  pressure  is  artificial  and 
denotes  the  nucleus  of  the  comet. 

Figures  3a  and  3b  show  the  plasma  density  in  the  comet 

magnetosphere.  The  perspective  in  the  sun-comet  -  IMF  (polar)  plane  is 
presented  in  3a,  and  3b  is  in  the  plane  sun-comet  perpendicular  Co  the  IMF 
(equatorial).  In  each  case  the  single  point  maximum  occurs  at  Che 
nucleus.  The  striking  feature  in  the  density  plots  is  Che  plasma  tail. 


The  tail  is  very  narrow,  ~  1  x  10^  km,  in  the  polar  direction  but  quite 

broad  in  the  equatorial  plane.  The  plasma  density  (assuming  CO^)  ranges 

—3  —3 

from  a  few  100  cm  near  the  nucleus  to  a  density  on  order  10  cm  at  the 

mid-line  of  the  tall  2  x  10^  km  behind  the  nucleus.  The  simulations 

indicate  that  the  plasma  density  structure  of  a  comet  is  a  spheroidal  head, 

a  few  times  10^  km  in  diameter,  with  a  long  ribbon-like  tail,  10^  km  thick 

in  the  direction  of  the  IMF  and  10^  km  wide  in  the  equatorial  plane.  This 

characteristic  plasma  structure  can  been  seen  more  clearly  in  cross- 

sectional  contours  on  the  plasma  tall  shown  in  Figure  4.  Figures  4a,  b,  c, 

and  d  display  plasma  density  contours  on  cross  sections  through  the  comet 

e  e  fL 

nucleus,  as  well  as  2.7  x  10  km,  8.3  x  10  km,  and  2  x  10  km  behind  the 
nucleus,  respectively.  The  cross-section  thru  the  nucleus  shows  circular 
contours.  As  the  cross-sectional  plane  is  moved  tallward  the  contours 
become  progressively  flattened  into  the  equatorial  plane.  At  distances 
further  into  the  tall  than  10^  km  the  contours  describe  a  rlbbon-llke 
plasma  structure  with  the  small  dimension  in  the  direction  of  the  IMF. 

In  the  preceedlng  results  we  have  shown  the  interesting  and  unique 
structure  of  a  comet  magnetosphere  and  plasma  tail  in  a  steady  solar 
wind.  The  results  are  essentially  in  a  steady  state  which  was  achieved 
after  simulating  13  hours  of  real  time.  The  stucture  of  the  comet  is 
especially  sensitive  to  the  solar  wind  IMF,  particularly  the  long  ribbon- 
like  plasma  tall.  In  the  remainder  of  this  section  we  will  present  results 
for  the  distortions  caused  in  this  structure  by  rotations  in  the  IMF. 

In  our  first  study  of  the  effect  of  a  non  steady  solar  wind  on  the 
cometary  magnetosphere  we  Induced  a  90°  rotation  of  the  IMF  in  the 
inflowing  solar  wind.  The  rotational  discontinuity  was  propagated  in  from 
the  front  of  the  mesh  as  a  one  cell  wide  rotation  of  the  IMF  from  the  z 
direction  into  the  y  direction.  The  rotational  discontinuity  required 
about  1  hour  to  propagate  to  the  cometary  bow  shock  but  changes  in  the 
cometary  tail  did  not  begin  to  show  clearly  until  the  disturbance  has 
propagated  far  past  the  nucleus  in  the  undisturbed  solar  wind.  Figures  5 
and  6  show  the  plasma  densities  6  hrs  and  figures  7  and  8  show  them  10  hrs 
after  introduction  of  the  rotational  discontinuity.  In  Figure  Sa,  the  old 
polar  plane  and  new  equatorial  plane,  we  see  the  central  remnant  of  the 
previously  existing  plasma  tall  and  the  formation  of  the  new  plasma  sheet 


in  the  outer  regions.  In  5b,  the  new  polar  plane  and  the  old  equatorial 
plane,  we  see  the  narrowing  of  the  plasma  tall  near  the  nucleus  and  the 
remnant  of  the  old  plasma  tail  far  from  the  nucleus.  Figures  6a  and  6b 
show  more  clearly  how  this  process  takes  place.  In  6a,  which  Is  closer  to 
the  nucleus,  the  plasma  tall  Is  twisting  Into  Its  new  orientation  from  the 
outside  In,  while  In  6b  further  from  the  nucleus  the  plasma  tall  still  has 
Its  original  orientation.  Figures  7  and  8,  10  hours  after  Introduction  of 
the  rotational  discontinuity,  show  a  continuation  of  the  process.  As  the 
rotational  discontinuity  propagates  Into  the  cometary  magnetosphere  at  the 
plasma  flow  velocity,  a  plasma  tall  In  the  new  orientation  gradually  Is 
created  and  propagates  Into  the  comet  behind  the  rotational  discontinuity. 
The  twist  In  the  plasma  sheet  shows  the  location  of  the  rotational 
discontinuity  as  It  moves  Into  the  cometary  magnetosphere.  The  flow  of 
plasma  In  the  comet  tall  remains  steady  and  In  the  anti-solar  direction 
throughout  the  twisting  of  the  tall  from  Its  original  position  to  Its  new 
orientation. 

In  our  second  study  of  the  effect  of  a  variable  solar  wind  on  the 
cometary  magnetosphere  we  Induced  a  180*^  field  reversal  of  the  IMF  In  the 
Inflowing  solar  wind.  The  results  for  -  the  field  reversal  are  shown  In 
Figures  9,  10,  and  11  for  times  of  6,  8,  and  10  hours  after  Introduction  of 
the  field  reversal,  respectively.  There  are  two  distinct  features  which 
appear  In  the  figures  of  the  plasma  density.  There  is  a  pronounced 
depletion  of  plasma  density  which  takes  place  at  about  10^  km  In  the  anti- 
solar  direction  from  nucleus.  The  maximum  depletion  is  about  a  factor  of  4 
and  the  depletion  propagates  tallward  with  the  local  plasma  flow 
velocity.  The  depletion  appears  to  arise  owing  to  small  changes  In  the 
plasma  flow  velocity  in  the  near  comet  tail  which  occur  shortly  after  the 
field  reversal  reaches  the  shock  front.  The  other  pronounced  feature  Is 
the  formation  of  ray-llke  structures  In  the  plasma  tall.  The  ray-llke 
structures  fold  with  the  magnetic  field  onto  the  tail  plasma  sheet  and 
reenhance  the  plasma  tall  densities  to  the  earlier  values.  The  ray 
structure  In  the  plasma  density  occurs  at  the  neutral  line  of  the  IMF 
reversal  and  may  be  associated  with  a  quiet  form  of  magnetic 
reconnection.  The  plasma  In  the  rays  flows  with  the  local  convection 
velocity  which  Is  a  l-sun'  .d  and  Is  virtually  undisturbed  from  the  steady 
values  existing  before  he  event. 


Discussion 


The  results  presented  above  provide  a  rather  simple  conceptual  model 
of  a  cometary  magnetosphere  and  Ion  tall.  The  solar  wind  carrying  an 
Imbedded  IMF  encounters  the  comet  as  a  massive  source  of  ions.  The 
cometary  Ions  decelerate  the  solar  wind  and  a  weak  shock  Is  formed  in  front 
of  the  comet.  The  IMF  field  lines  closest  to  the  source  of  Ionization 
"hang  up"  and  are  draped  over  the  comet  creating  a  long  magnetic  tall.  As 
this  process  occurs  the  solar  wind  has  barely  deviated  from  the  anti-solar 
direction  and  flows  essentially  straight  on  through.  On  each  IMF  field 
line,  the  region  nearest  the  source  of  Ionization  Is  retarded  the  most  and 
therefore  continues  to  accumulate  more  Ionization.  After  passage  through 
the  shock  and  cometary  head  the  plasma  reaccelerates  and  rarefies,  but 
those  regions  passing  closest  to  the  production  source  and  accumulating  the 
most  plasma  also  reaccelerate  and  ratify  more  slowly.  This  process  leads 
to  a  cometary  plasma  structure  with  a  more  or  less  spherical  head  and  a 
long  tenuous  rlbbon-llke  tall:  the  ribbon  tall  being  everywhere 
perpendicular  to  the  local  IMF.  Rotations  In  the  IMF  lead  to  twists  of  the 
rlbbon-llke  tall  with  the  inner  most  tail  regions  being  tied  to  the  head  by 
the  old  IMF  and  the  outer  regions  of  the  tall  oriented  to  the  new  IMF 
direction.  IMF  reversals  In  the  simulations  show  a  rarefaction  of  plasma 
In  the  tall  which  could  be  Interpreted  to  be  disconnection  events  and  also 
show  the  formation  of  ray-like  structures  which  fold  onto  and  reenhance  the 
tall  plasma  density.  We  will  now  discuss  certain  features  and  processes  In 
more  detail. 

The  simulation  results  strongly  argue  for  weak  magnetic  fields  In 
comet  tails.  There  has  been  considerable  disagreement  on  this  subject  In 
the  past  with  Hyder  et  al.  (1974),  Ip  and  Mendls  ( 1976a, b),  and  Mendls 
(1977)  arguing  for  strong  fields,  and  Ershkovlch  (1978,  1982)  and  the 
simulation  by  Schmidt  and  Wegmann  (1980,  1982)  arguing  for  smaller 
values.  Our  simulation  results  are  definitely  In  agreement  with  the  work 
of  Ershkovlch  as  the  maximum  field  In  the  head  Is  slightly  less  than  a 
factor  of  10  greater  and  at  distances  of  4  x  10^  km  In  the  tall  are  less 
than  a  factor  of  2  greater  than  the  IMF  value.  We  are  unaware  of  plasma 
kinetic  mechanisms  not  Included  In  the  MHD  formalism  which  could  be 
responsible  for  enhancing  the  magnetic  field  In  the  cometary  magnetosphere. 
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The  results  have  pointed  out  two  mechanisms  which  can  Lead  to 
formation  of  ray-llke  structures  In  the  comet  Ion  tall.  For  the  first 
mechanism  associated  with  a  reversal  In  the  IMF,  there  Is  an  accumulation 
of  plasma  In  the  comet  magnetosphere  along  the  IMF  neutral  line.  These 
plasma  accumulations  form  straight  ray-llke  structures  attached  to  the  head 
of  the  comet.  Whether  this  accumulation  Is  associated  with  magnetic 
merging  In  the  numerical  model  owing  to  numerical  resistivity  Is  at  present 
unclear.  We  hope  to  study  this  question  In  more  detail  In  the  future. 
Nevertheless,  the  plasma  enhancements  do  occur  In  the  results,  they  fold 
slowly  onto  the  tall  axis  and  enhance  the  plasma  density  there  forming  a 
straight  thin  plasma  tall. 

The  other  mechanism  Involves  the  helical  twisting  of  the  thin  rlbbon- 
llke  plasma  tall  by  rotational  discontinuities  In  the  IMF  of  less  than 
180*^.  An  observer  viewing  this  twisted  structure  from  the  side  would  see  a 
bright  ray  at  positions  where  the  rlbbon-llke  tall  was  aligned  with  the 
viewing  direction.  These  rays  created  by  viewing  perspective  would  also 
fold  onto  the  central  axis  appearing  to  brighten  the  tall  there  also.'  This 
mechanism  Is  similar  to  one  described  by  Schmidt  and  Wegmann  (1980,  1982), 
however  we  are  unable  to  ascertain  from  their  description  If  we  have 
Identified  the  same  mechanism. 

The  thin  ray-llke  structures  and  plasma  tall  which  we  Identify  in  the 
simulation  results  all  have  a  width  of  a  few  times  10^  km.  This  mlmimum 
size  Is  due  solely  to  the  resolution  In  the  numerical  grid.  Structures 
formed  In  a  higher  resolution  simulation  or  In  a  real  comet  could  be  much 
thinner . 

Another  feature  of  the  simulation  results  for  an  IMF  reversal  Is  the 
suggestion  of  a  plasma  tall  disconnection  event  as  proposed  and  discussed 
by  Nledner  and  Brandt  (1978,  1979)  and  Nledner  et  al.  (1981).  In  the 
simulation  results  we  observe  a  decrease  by  about  a  factor  of  4  In  the 
plasma  density  which  takes  place  at  an  anti-solar  distance  of  about  10^  km 
from  the  nucleus.  This  plasma  depletion  occurs  roughly  4  to  5  hours  after 
the  IMF  reversal  reaches  the  comet.  It  is  caused  by  small  flow  velocity 
changes  In  the  near  comet  plasma  sheet  which  are  possibly  associated  with 
pressure  changes  In  the  bow  region.  The  plasma  density  decrease  Is  not 
associated  with  magnetic  reconnection  In  the  plasma  tall.  Neither  Is  It 


directly  associated  with  the  IMF  reversal  region  which  is  responsible  for 
the  appearance  of  the  the  ray  structures  off  the  tail  axis.  The  depleted 
density  region  propagates  in  the  anti-solar  direction  at  the  local  plasma 
flow  velocity,  and  the  depleted  plasma  is  finally  replaced  considerably 
later  by  folding  of  the  rays  onto  the  tail  axis. 

Since  the  same  numerical  code  used  for  these  simulations  has  also  been 
used  for  simulations  of  the  earth's  magnetosphere  (Brecht  et  al.  1981, 
1982),  it  la  possible  to  make  some  interesting  comparisons  between  the  two 
sets  of  results.  In  the  case  of  the  earth  the  solar  wind  forms  a  strong 
shock  in  front  of  the  magnetosphere  and  deviates  to  flow  around  the 
magnetic  obstacle,  while  in  the  case  of  the  comet  the  shock  is  weak  and  the 
solar  wind  flows  essenc..'lly  straight  on  through.  In  the  case  of  the  earth 
the  results  show  a  distinct  magnetopause  with  a  clear  separation  of  solar 
wind  plasma  and  flow  from  that  of  the  earth's  plasma,  while  in  the  case  of 
the  comet  there  is  no  distinct  separation.  In  the  earth's  magnetotail  the 
plasma  is  extremely  tenuous  with  plasma  densities  less  than  1  cm  ,  while 
in  the  comet  tall  plasma  densities  are  larger.  The  plasma  sheet  of  the 
earth  separating  the  tall  lobes  is  a  tenuous  hot  plasma,  while  in  a  comet 
the  plasma  sheet  is  a  higher  density  cold  plasma.  In  the  case  of  the  earth 
disturbances  in  the  solar  wind  can  lead  to  dynamic  magnetic  reconnection 
with  high  velocity  plasma  flows  and  sudden  energetic  acceleration  events, 
while  in  the  case  of  the  comet-solar  wind  interaction,  disturbances  may 
lead  to  magnetic  reconnection  but  the  reconnection  is  slow  and  its  effects 
are  mild.  It  is  clear  that  the  two  magnetospheres  and  the  plasma  processes 
in  them  are  very  different. 

The  reader  might  wonder  at  the  remarkably  different  behavior  of 
magnetic  reconnection  in  Jhe  earth's  magnetosphere  and  that  of  a  comet.  At 
least  a  partial  answer  is  directly  related  to  the  orders  of  magnitude 
difference  in  the  Alfven  speed  in  the  two  plasma  tails.  In  the  earth's 
magnetosphere  the  Alfven  speed  can  approach  the  speed  of  light  near,  and 
possibly  even  in,  the  magnetic  neutral  sheet  region.  In  the  neutral  sheet 
region  of  the  comet  tails  the  Alfven  speed  is  orders  of  magnitude  lower. 
Since  the  speed  at  which  plasma  flows  out  of  the  reconnection  region  is 
controlled  by  the  Alfven  speed  in  the  inflow  region  (Petschek,  1964),  the 
plasma  flow  in  cometary  reconnection  events  is  bound  to  be  slow.  In  other 


words,  the  magnetic  field  energy  annihilated  by  reconnection  must  be  given 
to  the  plasma.  Since  the  8  of  the  cometary  plasma  is  >  1  to  begin  with, 
reconnection  cannot  lead  to  drastic  changes  in  the  plasma  behavior. 


In  conclusion  we  would  like  to  reiterate  the  main  features  of  the 
results.  The  solar  wind  interaction  with  a  comet  leads  to  a  magnetosphere 
and  magnetic  structure  as  suggested  by  Alfven  (1957).  The  magnetic  field 
in  comet  tails  is  weak  as  was  shown  by  Ershkovich  (1978,  1982).  The  plasma 
density  has  a  unique  structure  namely  a  spheroidal  head  stretched  out  in 
the  anti-solar  direction  into  a  long  ribbon-like  tail  Rotational 
discontinuities  in  the  IMF  lead  to  a  helical  twisting  of  the  ribbon-like 
tall  which  can  appear  as  rays  in  perspective.  The  results  for  the  IMF 
reversals  suggest  occurrence  of  tail  disconnections  as  described  by  Nledner 
and  Brandt  (1978,  1979)  and  lead  to  ray-like  structures  which  fold  onto  the 
tall  axis  to  form  a  new  plasma  tall. 
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Fig.  1.  Contours  of  magnetic  field  strength  and  vector  direction  in  the  sun-comet 
IMF  plane.  Contour  levels  of  field  strength  are  2-7.6^;  3-11^,  4  -  16,, 
S  -  22,,  6  •  29,,  7  -i  38,,  and  8  -  48,.  The  comet  is  located  at  (0,0)  and  dis¬ 
tances  along  the  axis  are  in  km. 
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Fig.  2.  Plasma  flow  velocity  vectors,  (a)  and  perspective  plot  of  plasma  pressure, 
(b).  The  flow  velocity  in  the  solar  wind  at  the  left  edge  is  400  km  sec~'  and  in  the 
central  comet  is  about  10  km  sec'*.  The  flow  shows  little  deviation  around  comet 
and  is  essentially  straight  on  through.  The  left  scale  on  the  pressure  plot  is  the  log* 
arithm  to  the  base  10  of  pressure  in  eV  cm~^.  On  the  pressure  plot  note  the  large 
pressure  jump  in  front  of  and  at  the  shock;  also  note  the  large  pressures  in  the 
magnetic  fleld  null  region.  The  single  point  maximum  is  the  location  of  the  comet 
nucleus  and  is  artificial.  Both  plots  are  in  the  sun-comet,  IMF  plane. 
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Fig.  4.  Plasma  density  contours  on  cros8>sectional  planes  through  the  comet  and  tail 
perpendicular  to  the  solar  wind  flow.  Diagram  (a)  shows  a  plane  through  comet 
nucleus,  (b)  is  2.7  x  10*  km  in  the  anti  solar  direction  from  the  nucleus,  (c)  is 
8.3  X  10*  km  in  the  and  solar  direction,  and  (d)  is  2.  x  10*  km  in  the  and  solar 
direcdon.  Notice  that  the  plasma  density  contours  are  nearly  circular  in  the  comet 
head  but  become  progressivehr  flattened  onto  the  equatorial  plane  further  into  the 
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Fig.  4  (Continued  Plume  density  contours  on  cross-sectionel  planes  through 
the  comet  and  tail  perpendicular  to  the  solar  wind  flow.  Diagram  (a)  shows  a 
plane  through  comet  nucleus,  (b)  is  2.7  x  10*  km  in  the  anti  solar  direction 
from  the  nucleus,  (c)  is  8.3  x  10*  km  in  the  anti  solar  direction,  and  (d)  is 
2.  X  10*  km  in  the  anti  solar  direction.  Notice  that  the  plasma  density  contours 
are  nearly  circular  in  the  comet  head  but  become  progressively  flattened  onto 
the  equatorial  plane  further  into  the  tail. 


(b) 


Fif.  6.  Plasma  density  contours  in  cross-sectional  planes  perpendicular  to  the 
solar  wind  6  hours  after  the  90*  routional  discontinuity.  Figure  (a)  shows  a 
oosa-sectional  plane  8.3  x  10^  km  behind  the  comet  nucleus  and  (b)  shows  a 
cross-sectional  2  x  10*  km  behind.  Note  the  twisting  of  the  plasma  sheet  from 
outside  inward. 
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